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Abstract 

The RP 2 gauge model which allows interpolation between the RP 2 and 0(3) spin 
models is studied in 2D. We use Monte-Carlo renormalization techniques for blocking 
the mean spin-spin interaction, {A}, and the mean gauge field plaquette, (P). The 
presence of the 0(3) renormalized trajectory is verified and is consistent with the known 
three-loop /3-function. The first-order 'vorticity' transition observed by Solomon et al. 
HI is confirmed, and the location of the terminating critical point is established. New 
scaling flows in ({A), (P)) are observed associated with a large exponent k in the range 
4 ~ 5. The scaling flows are found to give rise to a strong cross-over effect between 
regions of high and low vorticity and are likely to induce an apparent signal for scaling 
in the cross-over region which we propose explains the scaling observed for RP 2 and 
RP 3 models and also in a study of the SO (4) matrix model [||. We show that the 
signal for this 'pseudo' scaling will occur for the RP 2 spin model in the cross-over region 
which is precisely the region in which computer simulations are done. We find that the 
RP 2 spin model is in the same universality class as the 0(3) spin model but that it 
is likely to require a very large correlation length before the true scaling of this class 
sets in. We conjecture that the scaling flows are due either to the influence of a nearby 
new renormalized trajectory or to the ghost of the Kosterlitz-Thouless trajectory in the 
associated XY model. In the former case it is argued that the 'vorticity' fixed point 
controlling the critical behaviour terminating the first-order line cannot be identified 
with the conjectured new renormalized trajectory. 
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1 Introduction 



The nature of the phase diagram for two-dimensional RP models has been the 
subject of much recent discussion ||, 0. In ||, Caracciolo et al. compare the 
correlation length computed from simulation with that predicted from the per- 
turbative (3— function using the exact results for the mass gap in O(N) models. 
They found that for RP 2 (RP 3 ) the observed correlation length on lattices up 
to L = 512 was smaller than the expected value by a factor of 10 7 (10 4 ). Their 
conclusion was that either the asymptotic regime is very far indeed removed from 
the regime of their study, requiring lattices of sizes of 10 9 (10 5 ), or that these 
theories were not asymptotically free but that there exists a phase transition at 
finite f3 (non-zero temperature). Caracciolo et al. indeed provide evidence for 
the latter scenario by showing that their data scale in a manner consistent with 
a Kosterlitz-Thouless parametrization. The two persuasive features are thus that 
the correlation length is much smaller than that expected assuming an asymp- 
totically free theory and that scaling of the data is observed. This phenomenon 
occurs in a large class of models and the question is whether the signal for a phase 
transition at finite f3 and the observed scaling of data are genuine or not. 

The same effects have been observed to a less extreme extent by Hasenbusch 
and Horgan |J who investigated the continuum limit of the SO (4) matrix model. 
The measured ratio of the mass gap to A.^jg was compared with the theoretical 
prediction obtained using the Bethe ansatz ||. There was a disagreement be- 
tween theory and experiment by about a factor of four the measured correlation 
length being about four times smaller than expected. However, the measurement 
using the covering group was in excellent agreement with theory. The numerical 
method used, due to Luscher et al. 0, relies in part on measuring the correlation 
length in a large volume and establishing that scaling holds with only small and 
perturbative violations. Although in the 50(4) case there were strong indica- 
tions that the results scaled, the discrepancy between simulation and theory led 
to the conclusion that the signal for scaling was only apparent and that a true 
continuum limit had not been achieved in the large volume simulation. It was 
conjectured that the cause of the deception was the presence of vortices in the 
SO (4) model which are absent in the case of the covering group since 

n 1 (50(4)) = Z 2 , ^(5^(2)) = . (1.1) 

One question is, therefore, whether a bogus signal for scaling can be observed 
in the presence of vortices in two dimensions. In the work presented here this 
question is addressed in the context of an RP 2 gauge theory which allows an 
interpolation between the pure RP 2 and 0(3) spin models. This gauge model 
contains Z 2 vortices coupled to a chemical potential. We observe the conventional 
0(3) renormalized trajectory and show that our results are consistent with the 
known three-loop f3 function. We establish the existence of a first-order transition, 
first suggested by Solomon et al. JTJ , for which the order parameter is the vorticity. 
The critical point terminating this first-order line will be in the domain of a 
new 'vorticity' fixed point. Using Monte-Carlo renormalization group (MCRG) 
techniques we observe certain flows on which the blocked observables scale and 
suggest that these scaling flows are due to the influence of a nearby renormalized 
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trajectory which gives rise to the possibility of the existence of a fixed point 
other than the 0(3) one. We argue that it is unlikely that any new fixed point 
can be identified with the inferred 'vorticity' fixed point. Our results strongly 
indicate that the apparent or 'pseudo' scaling behaviour is due to a cross-over 
effect associated with the proximity of the new scaling flows to the line of RP 2 
spin models in coupling constant space. The cross-over is between regions of high 
and low vorticity which emphasizes the crucial role of vorticity in the observed 
properties of the model. Where relevant our results confirm or complement those 
obtained by Solomon et al. in an earlier study of this model. 

Another reason for studying the RP 2 gauge models is that it has been conjec- 
tured || that in 2D the continuum limit in the RP 2 spin model is distinct from 
that in the 0(3) spin model. Niedermayer et al. @ and Hasenbusch |§ have 
suggested that this conjecture is incorrect and that there does exist a continuum 
limit in the RP 2 model which is controlled by the 0(3) fixed point. The essential 
question is whether or not the RP 2 model is in the same universality class as the 
0(3) model. By using MCRG methods to show the topology of renormalization 
group trajectories in the RP 2 gauge theory we find that a consistent and simple 
interpretation of our results is that the RP 2 and 0(3) models are in the same uni- 
versality class: an interpretation which supports the conclusions of Niedermayer 
et al. [0] and Hasenbusch ||. 

All results are for RP 2 gauge models but the simulation can be generalized to 
jlpN-i anc j a cursor y investigation for N > 3 has indicated that broadly similar 
results hold for this general case. 

In section |2] we define the model under study; in section |3| we briefly describe 
the simulation techniques; in section |] we define the Monte-Carlo renormalization 
group method used and describe the measurement procedure; in section [5] we 
present the results and discussion and in section [7| we draw our conclusions. 



2 The model 

The action used is 

5({5},{s}) = -p[Y,S x -S x+fl a x ^ + h'EPm) , 

\£E,/i X / 

where x = (xi,X2) , x±,X2 G Z, 1 < x±,X2 < L, labels the sites of an integer 2D 
square lattice of side L, and takes values in = (0, 1), /x 2 = (1; 0) • The spin 
S x is a unit length three-component vector at site x and <r X: M is a gauge field on 
the link (x, /x) taking values in [1, —1] . The plaquette of gauge fields is denoted 
by P x {<y) where 

Px(@} ^x,/j, 1 (^x+fi 1 , /x 2 &x+(j, 2 , Mi < ^" a; ! A*2 ' 

This action is invariant under the gauge transformation 

&x * 9x & x , 
a x,^ ~^ 9x a x, fj, 9x+n i (2-2) 
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with g x e [1, -1] . 

Vortices reside on plaquettes where P x (a) = — 1, and are suppressed (en- 
hanced) if the chemical potential, /z, is positive (negative). The pure 0(3) model 
corresponds to fi — > oo and the pure RP 2 model corresponds to fi = . 



3 The simulation 

A local update was used comprising a combination of heat-bath, microcanonical 
and demon schemes. For fixed gauge fields the spins {S} were first updated by 
a heat-bath algorithm which can be generalized to 0(N) and so for this section 
we will consider S x to be an N— component spin of unit length. The heat-bath 
method is to project each spin onto a 3D subspace of the N- dimensional space 
in which the spins take their values. The 3D subspace is chosen at random but 
is the same for all spins during one lattice update. Let the projection of S x onto 
this space be denoted R x . Then clearly 

{R x )i = (S x ) n , 2 = 1,2,3, 1 < jt < j 2 < j 3 < N , (3.1) 

where the ji are chosen randomly subject to the restrictions above. The single-site 
probability distribution for R x is then 

Q(R X ) ex exp(M x ■ R x ) , (3.2) 

where 

Vx-^ti) ■ ( 3 - 3 ) 

M 

The heat-bath update at each site is done by replacing R x by R' x which is 
chosen with the probability distribution Q(R X ). The microcanonical spin update 
is 

2 (S x ■ M x ) M x . , 

S x - S' x = - S x + 1 * g x - . (3.4) 

I x I 

The demon update is applied to the gauge fields only. In general, it is only 
necessary to introduce one demon variable for the whole gauge configuration. 
However, when running on a massively parallel computer it is necessary to have 
one demon per processor and then each demon must migrate through the whole 
lattice. This is easily achieved by moving demons sequentially between processors. 
We illustrate the method with one demon variable d, d > 0. The action in 
equation (^) is augmented by the demon to become 

S demon ({S},{a},d) = S({S},{a}) + (3d. (3.5) 

Then for each link (x, ii) the trial gauge field update is (cr x ^ , d) — > (—cr X fl , d') 
where d! is chosen so that Sdemon is unchanged. That is, 

d' = 

d 2 <J X ^ { S x • S'cc+zx fJ> i®x, i> @x+i>, [A @x+n, i> @ x— v, v (fx—v, n (^x—v+iJ,, u) } ) 

(3.6) 
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where v is the orthogonal vector to /x. The update is accepted only if d' > 0. 
Note that the update is microcanonical in the augmented configuration space of 
fields plus demon and hence it is independent of (3. 

One complete lattice update consisted of one heat-bath update followed by 
an alternating sequence of Nmd microcanonical and demon updates. The value 
of Nmd that optimizes the decorrelation of the configurations depends on many 
factors and we did not spent much effort in tuning Nmd but regard Nmd ~ 10 
as a reasonable value. The heat-bath update took about ten times the time of 
the combined microcanonical and demon updates and so there was little time 
penalty for this choice. Depending on the coupling constant values we found that 
decorrelated configurations were produced within 2-30 iterations. Lattice sizes 
ranged from 64 2 to 512 2 , and typically the numbers of configurations per run 
were e.g., 2 • 10 6 for 64 2 and 5 • 10 5 for 256 2 . 

The simulations were carried out on the HITACHI SR2201 computers in the 
Cambridge High Performance Computing Facility and in the Tokyo Computing 
Centre. 



4 The Monte- Carlo renormalization scheme 

The objective is to establish the topology of renormalization group flows in the 
relevant large-scale variables and infer the phase structure of the model. After 
sufficient blocking we assume that we are dealing with renormalized observables, 
and so different phases will be distinguished by singularities in the renormalization 
group flows. This has been discussed, for example, by Nienhuis and Nauenberg 
and by Hasenfratz and Hasenfratz We assume that there are at most 



two relevant couplings in the neighbourhood of any fixed point in which we are 
interested. We also assume that the chosen blocked operators have components 
which span the two-dimensional space of relevant operators, i.e., the operators 
conjugate to these relevant couplings. From our earlier experience || and from 
the surmize stated in the introduction that vorticity plays a vital role, we chose to 
study how the mean values of the spin-spin interaction, A, and of the plaquette, 
P, flow under blocking. For a given configuration these quantities are defined by 
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A - E S »° ' s * 



W x,, 



X+fJb (I ) 



p = ~ E p *(°) ■ (4- 1 ) 

* X 

(A) lies in [0,1] and (P) lies in [—1,1] and the mean vorticity is defined by 
V = (l-P)/2. 

For each configuration {S, a} on a lattice of side L we derive a blocked con- 
figuration {S B ,a B } on a lattice of side L/2 . The blocking transformation for 
the spins is 

S x B = 

I numerator I 
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(4.2) 



Based on earlier work by Hasenbusch et al. [ I Ij the parameter a was chosen to 
be 0.0625 . Choosing other reasonable values for a was found not to change any 
outcome or conclusion. This gauge-invariant blocking transformation is shown in 
figure |l| 

To block the gauge field the products of gauge fields were computed for the 
three Wilson lines joining the end points of the blocked link shown in figure [l]. 
These field products are denoted by Wq, W + , W_. For the blocked link joining x 
to x + 2/j, the Wi are given by 

&X, fl&X + fl, fl ) 

W+ 

— @x, u Gx+v, fj, ^X+fJL+U, fj, &X+2fJ,, U 1 

W— — °~x-v, v Vx-v, /i, a x+n-h>, n °~x+2fi-is, v > (4.3) 

where v is the orthogonal vector to \x. The blocked gauge field was assigned the 
majority sign of the Wf 

B = w + + w + w„ 

*b,hb \W+ + W + W-\ ' V ' ; 

This blocking transformation has the important property that it ensures that two 
vortices on adjacent plaquettes of the original lattice will cancel and not survive in 
the blocked lattice. This is clearly true if the vortices lie in the same 2x2 block 
since they add mod2, but the majority rule guarantees cancellation also when 
two adjacent vortices lie either side of the block link separating two neighbouring 
blocks. This is illustrated in figure |[ 

For a given pair of coupling constants (/3, fi) and given lattice size L x L, 
each configuration was blocked by successive transformations until the blocked 
lattice size was 8x8. The operator expectations (A) L (f3,ii) and (P) l((3, /z) were 
then measured and averaged over all configurations. For given (j3, fi) this was 
done for L = 64, 128, 256, 512 which gives four points on a segment of a flow in 
the ((^4), (P)) plane with each point labelled by initial lattice size. Each point 
corresponds to a rescaling of length by a factor of two compared with the previous 
point. 

The errors on the observables were determined by averaging the results for 
successive configurations in bins of 2 M , M = 0,1,2,... and calculating the er- 



rors on the ensemble of bin-averaged measurements [JT^]. The true error is the 
asymptotic value achieved for large enough M. The decorrelation length can also 
be estimated from the behaviour of the error as a function of M. The number 
of independent configurations ranged from about 600 for L = 512 to in excess of 
2 • 10 4 for L = 64 . Errors were also estimated from the ensemble of independent 
measurements from different processors. 

5 Results 

Each flow segment consisting of four points, but more complete flows can be built 
up by extending the flow in either direction by tuning to new couplings (/?', //) 
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so that 



{P)v{?,lt) = (P)l(P,I*), (5-1) 

for some L and V . The flow for given (/?',//) can then be computed. In general 
this will be an approximate procedure because a segment in the ((A), (P)) plane is 
the projection onto this plane of part of a full flow in the higher dimensional space 
of observables. By tuning as described we are ensuring only that the projections 
of blocked points coincide not the blocked points themselves. In principle, we 
need to match a full complement of observables by tuning a complete set of 
couplings, conjugate to these observables, which define the most general action 
consistent with the symmetry. In general, the effect of couplings which are not 
included cannot be properly taken into account. However, we assume that in the 
neighbourhood of a fixed point there will be at most two relevant couplings and 
that the projection onto the ((3, //) plane of the space they span is non-singular. 
The effect of irrelevant couplings is mitigated by performing an initial blocking 
by a factor which will significantly reduce the errors induced by the projection 
so that we are effectively dealing with renormalized operators. Since the target 
lattice is always 8x8 the size of this factor depends on the initial lattice size and 
hence on available CPU. For our study this initial blocking factor had a minimum 
value of eight. 

There are errors due to finite-size effects which can be parametrized in terms 
of the parameter z = £l/L where £x, is the correlation length on the Lx L lattice. 
Because the target lattice is the same size throughout, the values of z associated 
with coinciding points, equation ( |5.1|) , are similar and so the mismatch in the 
finite-size errors between different segments joining up to make a longer flow will 
be minimized. There will nevertheless be a residual finite-size effect which it is 
generally not possible to estimate except in the case of the 0(3) spin model which 
is discussed in the next section. 

The details of the flows can depend on z and the details of the blocking 
scheme. It follows that conclusions about the physical properties of the theory 
can be deduced only from universal or topological properties of the flows such as 
the occurrence of fixed points and singular behaviour. 



5.1 The 0(3) renormalized trajectory 

In the limit /i^oowe recover the 0(3) spin model, and to test our procedures 
and assumptions we should, at the very least, be able to recover the perturbative 
/^-function for this model. The projection of the 0(3) renormalized trajectory 
onto the (/3, /i) plane is the /3-axis. We expect corrections to scaling due to finite 
lattice-spacing artifacts which will be a function of L. We find that our method 
works well for sufficiently large j3 once the tree-level approximation for these 
scaling corrections has been taken into account. Consider the block observable 

1 ~ ' °y B 

large (3 we write 



A = S Xg ■ S yB with the blocking transformation defined in equation |4.2| . For 



S x = qJl-(3-i<f>l + (3~ 1/2 <t> x , (5.1.1) 
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where q = (1, 0, 0) and <fi = (0, fa, fa). Using [O] and keeping terms up to j3 1 
we find after one blocking step that 



xb 



with 



+ </w M 



(5.1.2) 
(5.1.3) 



XB ' l + Aa 

and where \x is summed over nearest neighbour links. For large f3 we expect the 
blocked expectation value, (A)l, on the 8x8 lattice to behave as 

(A) L = 1 - C(L)/P + 0(l/(3 2 ) . (5.1.4) 

For large enough L we expect C(L) to attain its limiting value. However, there 
is still some variation in C(L) for the values of L we are using. In order to 
accommodate the bulk of this correction to scaling we define the effective coupling 
Peff by 

and modify the matching condition of equation ( |5.1| ) in this case to become 

P eff (L,P) = f3 eff (L>,P>) . (5.1.6) 

We then expect that 



log(L/L') 



du 



(5.1.7) 



where u = 1//3, u' = 1/(3'. 

C (L) is determined from a free field theory calculation on an L 2 lattice of the 
kinetic term for the blocked field <& x which is defined on the target L 2 B lattice by 
iteration of equation ( |5.1.3|) . This calculation is done easily numerically and the 
results are given in table [I]. We use L B = 8 in subsequent calculations. 



L 


C(L,L B = 8) 


C(L,L B = A) 


16 


0.53890071 


0.54363018 


32 


0.56795927 


0.55935100 


64 


0.58370587 


0.56866918 


128 


0.59303041 


0.57453025 


256 


0.59889307 


0.57829938 


512 


0.60266260 


0.58074276 



Table 1: The function C(L) defined in equation ( 5,1.4| ) for blocking from an L 2 lat 



tice for target lattices with L B = 8,4 . C(L) gives the tree approximation for the 
dependence on L/a of corrections to scaling in the 0(3) spin model. 

Using equation ( |5.1.7|) and the three-loop /3-function from |J we determine 
sequences for the bare coupling, /3, for which successive terms correspond to 
blocking by a factor of two. In table || we compare the values of /3 e ff for the 



8 



two sequences (3 = 5.0, 4.8861, 4.7721 and 2.0, 1.8803, 1.7560. If our simulation 
reproduces the correct /3-function then the matching condition 

P eff (2^L,p n ) = /3 eff (L,p n ), (5.1.8) 

must be satisfied, where (3 n is the n-th term in the sequence. From table |2| we 
see that this condition is indeed very well satisfied for the sequence starting with 
(3 = 5.0. For the other sequence the larger values of /3 e ff agree well and only 
as (3 e ff decreases is there an increasing discrepancy which signals a significant 
deviation of the three loop approximation to the /3-function from the correct 
value and also the possible effect of the neglected L-dependent 0{1/ (3 2 ) terms 
in equation ( |5.1.5|) . Our expectation is confirmed that the method correctly 
reproduces asymptotic scaling and probes the renormalization group flow close 
to the renormalized trajectory. 



p 


Initial lattice size, L 


64 


128 


256 


512 


5.0 




4.355(1) 


4.227(3) 


4.107(5) 


4.8861 


4.343(1) 


4.244(3) 


4.110(3) 


3.995(5) 


4.7721 


4.249(1) 


4.124(1) 


4.002(3) 




2.0 




1.2752(4) 


1.1366(4) 


0.9852(7) 


1.8803 


1.2733(2) 


1.1289(4) 


0.9729(7) 


0.8194(4) 


1.7560 


1.1210(3) 


0.9532(4) 


0.8000(2) 





Table 2: Values of f3 e ff(L, f3) for sequences of bare coupling j3 computed using equation 
( 5.1.7 ) where successive couplings in the sequence correspond to blocking by a factor of 
two. The required matching condition, equation ( 5.1. 8| ) , for verification of asymptotic 
scaling is very well satisfied for the sequence at larger j3 and the deviation in the other 
sequence is largest for the smaller values of (3 and is due to deviation of the three-loop 
approximation to the /3-function from the true value and possibly to terms neglected 
in equation ( [5.1.5 ) . 



5.2 New scaling flows 

In figures |^ to H we plot the flow segments for various (/?, /x) values in the ((A), (P)) 
plane, where the longer flows in figures |3| and [| are composed of superimposing 
segments using equation fl5.1| ). There is a flow on which the observables scale. 
This is shown in figures |3] and |4j. Nearby flows also showed scaling but are not 
included in the figures for reasons of clarity. To see that observables scale each 
flow of four points was successively overlaid using the tuning described in equation 
( |5.1|) with L' = L/2. For the scaling flows the points of the overlaid flow segments 
coincide very well within errors (the tuning is not absolutely exact) showing that 
there are only very small effects from the irrelevant operators. In table we give 
the values of ((A), (P)) for points on the segments making up two such scaling 
flows, labelled set 1 and set 2, together with the initial lattice sizes and coupling 
constant values. 
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The tuning described by equation ( |5.1| ) was done by trial and error. We 
tried to estimate the effect of small increments in the coupling constant values by 
using the method of reweighting but this was found not to work. This was mainly 
because the blocked ((A), (P)) values were very sensitive to the initial couplings 
and so to use reweighting was not a realistic possibility. Also it was found that the 
effect of a change in (3 could be partly compensated by a change in \i. This was 
because the Jacobian d((A), (P))/d((3,n) was relatively small and presumably a 
better choice for the pair of observables and/or initial couplings would increase 
its value. However, in no case was this Jacobian dangerously small and tuning 
by trial was effective. 

For these flows we did not apply a compensation for scaling corrections of 
the kind used in the previous section for the 0(3) model. It is not possible to 
carry out a similar tree-level calculation to determine the compensation, if any, 
as a function of L/a since the model and the appropriate coupling(s) associated 
with the scaling flows are not known. Also, the data presented in this section 
were obtained before the details of the 0(3) analysis were known and the possible 
residual dependency on L analyzed. Overall, the deviation from scaling for the 
scaling flows is not large and we must allow for the possibility that it is due, in 
some measure, to the finite lattice-spacing artifact of the kind analyzed for 0(3). 
However, there are a few points shown which do deviate substantially from scaling 
for which such an explanation is unlikely. For both sets these points are for the 
two largest values used for f3, and for the largest initial lattice size, L = 512. We 
defer the discussion of the reasons why they do not scale well until section |TJ. 

For the couplings which generate the scaling flow denoted by set 1, a plot of 
exp(— /j,) versus g = 1/(3 can be approximated by a straight line with all except 
the one with smallest (3 well fitted by 

exp(-/i) = 1.74(1) - 2.98(4)3 , (5.2.1) 

Phenomenologically, we find that it is consistent to associate a scaling exponent, 
k, with the scaling flow using the relation 

k = plog(2)/log(( gi - g*)/(g 2 -g*)) . (5.2.2) 

where the number of blockings for g\ and g 2 differ by p (a factor of TP in change of 
scale) in order for the corresponding points in the ((A, ), (P, )) plane to coincide 
on the scaling flow. The value of g* is a free parameter which is chosen to obtain 
the best fit assuming k in equation ( 5.2.2| ) is constant. Even so, g* is not very 



well determined by this alone and we find a range of values g* ~ .13 ~ .15 to be 
acceptable, for which the various pairings of couplings for set 1 from table ^ give 
the values for k shown in table [|. The values of fi* in table |] are inferred using 
the linear fit of equation ( 5.2.1|) . Note that the possible values of /i* in table |] are 



far removed from the value \x = which characterizes the RP 2 model. Taking k 
in the range 4 ~ 5, the dimension of the relevant operator associated with the 
scaling flows is 

A = D-I/k « 1.75 ~ 1.8. (5.2.3) 
A word of caution. The definition of k used in equation ( |5.2.2| ) is appropriate 



for conventional second-order behaviour. However, the consistency of the fit 
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(/?,//) 


Initial lattice size, L 


64 


128 


256 


512 


Set 1 


4.26, -0.045 


0.81894(3) 
0.98330(4) 


0.80587(9) 
0.9809(2) 


0.7913(2) 
0.9753(3) 


0.7761(3) 
0.9684(8) 


4.0, 0.0 


0.80793(3) 
0.98021(6) 


0.7913(2) 
0.9740(3) 


0.7690(2) 
0.9561(4) 


0.739(1) 
0.927(3) 


3.72, 0.06 


0.79261(5) 
0.97402(9) 


0.7678(1) 
0.9547(2) 


0.7273(4) 
0.9053(8) 


0.6424(8) 
0.770(3) 


3.45, 0.125 


0.76947(6) 
0.9576(1) 


0.7251(2) 
0.9025(4) 


0.6363(4) 
0.7601(10) 


0.4995(6) 
0.558(1) 


3.18, 0.187 


0.7233(2) 
0.9025(4) 


0.6298(2) 
0.7528(4) 


0.4872(2) 
0.5430(4) 


0.2878(5) 
0.3348(5) 


Set 2 


4.00, 0.02 


0.81574(3) 
0.98787(5) 


0.80228(6) 
0.9849(1) 


0.7871(2) 
0.9782(4) 


0.7694(9) 
0.967(2) 


3.73, 0.0833 


0.80236(4) 
0.98518(6) 


0.7846(1) 
0.9768(2) 


0.7611(2) 
0.9571(5) 


0.720(1) 
0.905(3) 


3.495, 0.14 


0.78475(7) 
0.9769(1) 


0.7577(1) 
0.9538(3) 


0.7085(3) 
0.8865(7) 


0.611(2) 
0.727(3) 


3.3, 0.18 


0.7596(1) 
0.9544(2) 


0.7070(1) 
0.8826(3) 


0.6062(4) 
0.7178(8) 


0.4564(7) 
0.506(1) 


3.06, 0.23 


0.7050(3) 
0.8821(7) 


0.6002(30 
0.7105(7) 


0.448(5) 
0.4990(9) 


0.2383(2) 
0.2973(4) 



Table 3: Values of the blocked observables ((^4), (P)) for the flow segments constituting 
two neighbouring scaling flows labelled by set 1 and set 2. The values of the couplings 
/i) labelling each segment is given and the different points on a given segment are 
labelled by the initial lattice size L. The blocked observables were measured after 
blocking to a fixed target lattice size of 8 x 8. Scaling of the blocked observables can 
be seen to hold extremely well for each set by noting that the ((^4), (P)) values lying 
on any given diagonal sloping from bottom left to top right agree very closely indeed, 
except for the largest two values of f3 for L = 512. Set 1 is shown in figures |3| and |j. 



should be taken only as a phenomenological parametrization and it is likely that 
a different scaling form such as derived from Kosterlitz-Thouless behaviour would 
give an equally good fit. For example, in |L3j Seiler et al. study the Z(10) model 
in 2D which they presume has a KT transition. They find that both second-order 
and KT scaling forms give equally good fits near the transition and that, indeed, 
the KT form is hard to reconcile with conventional theory. 



5.3 The first-order transition 

Solomon et al. [0 use simple arguments to suggest that a line of first-order 
transitions occur in the range fii < /i(/5) < fi2, Hi = — 0.293, /i 2 ~ —0.26. This 
line of transitions will terminate in a critical point associated with a continuous 
transition implying that a critical surface intersects the (/3, /i) coupling constant 
plane. Confirmation of the first-order transition is shown in figures [| and |6] where, 
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(<?*, A 4 *) 


Scaling exponent k 


0.13,-0.305 


5.10 


4.92 


4.919 


4.90 


4.76 


4.83 


4.84 


4.91 


4.88 


4.84 


0.15,-0.261 


4.19 


4.10 


4.15 


4.18 


4.02 


4.13 


4.18 


4.25 


4.26 


4.28 



Table 4: The scaling exponent k calculated using equation ( 5.2.2| ) and various pairings 
of couplings g = 1/(3 from set 1 given in table ||. Two different choices for (g* , //*) are 
used which correspond to the range giving consistent results for k. 



respectively, the values of (A) and (P) for a lattice of L = 128 are shown plotted 
against [i for fixed /3 = 6.0, 7.0, 7.5, 8.0. It is clear that there is no transition for 
(3 = 6.0, 7.0 but that there is likely to be a first-order transition for (3 = 7.5, 8.0. 
This places the critical point in the range 7.0 < (3* < 7.5 which is consistent 
with the investigation of Solomon et al. [Q. The value of fi* varies a little but is 
close to —0.26. There is no detectable dependence on L and studies on lattices 
with L = 256, 512 have shown identical results within statistical errors. From 
figures [| we infer that a good order-parameter distinguishing the two phases is 
the vorticity. There is also a discontinuity in (A) shown in figure |B| which is, 
however, not independent of the discontinuity in (P). The value of (A) will vary 
rapidly since it is sensitive to, and thus reflects, the discontinuous change in the 
vorticity. The effective potential will show no discontinuity and it is crude but 
reasonable to suppose that one particular linear combination of A and P plays 
this role. In figure |7j we plot (P) against (A) and indeed it can be seen that 
there is no sign of a sharp discontinuity and that the locus of points is reasonably 
linear. The outcome is that the combined operator 

U c (x) = E^'^+m^/x + iPM , (5.3.1) 

with 7 —0.29, has a continuous expectation value across the transition. The 
orthogonal combination, Ud{x), has a discontinuous expectation value which is 
sensitive to the vorticity and is a good order-parameter. The simple and persua- 
sive argument of Solomon et al. |L[] is based on minimizing the energy at (3 = oo 
to show that a first-order transition occurs in (P) at //(oo) = l/y/2- I. The 
argument also depends on continuity of the energy which, at (3 = oo, means that 
Uc{x) is identified with the local energy operator. This will be only approximate 
for (3 < oo. 

The critical point terminating the first-order line will be in the domain of a 
fixed point, the 'vorticity' fixed point, with a renormalized trajectory on which 
the vortex density scales thus defining a new continuum theory with non-zero 
vortex density. 



5.4 Cross-over of flows 

In figures |8| and |9] there are two further sets of flow segments, sets 3 and 4, which 
each show a clear cross-over as a function of initial couplings. Set 3 lies to the 
left of set 4. These two sets are examples of the cross-over effect which we infer 
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occurs in a narrow region formed by the neighbourhood of a continuous line of 
theories in the (f3, /i) plane. 

The couplings associated with sets 3 and 4 are listed in table |5] where, for each 
set, the couplings reading from left to right label the segments in order from the 
lowest in the figure (lowest (P)) to the highest (highest (P))- 





(P, I*) 


Set 3 


3.3 


3.3 


3.3 


3.3 


3.3 


3.3 


0.34 


0.36 


0.40 


0.45 


0.50 


0.60 


Set 4 


3.9 


4.05 


4.17 


4.19 


4.29 


4.5 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



Table 5: The couplings associated with sets 3 and 4. For each set, the couplings reading 
from left to right label the segments in order from the lowest in the figure (lowest (P)) 
to the highest (highest {P)). 

We concentrate in particular on the flow segments of set 4 which correspond 
to the RP 2 spin model (/i = 0). These segments are shown for various (3 in 
the range 3.9 — 4.5 where the cross-over in the flows is very strongly marked, 
occurring between (3 = 4.17 and (3 = 4.19. The values of ((A,), (P, )) blocked 
from L = 64 — > L = 8 show little variation in values but when blocked from larger 
L the variation is very strong indeed. Moreover, it is clear that for (3 = 3.9 ~ 4.1 
the flow is dominated by the proximity of the scaling flows, represented as the 
dashed line in the figure, which will induce a signal for scaling in pure RP 2 when 
(3 ~ 4. However, as (3 is increased there is a strong crossover effect, and for 
(3 > 4.25 the flow is consistent with dominance by the renormalized trajectory 
associated with the asymptotically free 0(3) fixed point at (f3*, n*) = (oo,oo). 
Consequently, the apparent scaling signal will only be transitory. For finite (3 
there will always be some free vortices, but for (3 > f3 cross - over we expect that the 
density of free vortices will vanish faster than l/£ 2 and the vortex density will 
not scale. This behaviour is consistent with the absence of a phase transition at 
finite (3 as well as with the 0(3) type continuum limit in RP 2 . 



6 Discussion 

The observation of the scaling flows reported in section |5.2| poses the question of 
whether we can attribute them to the influence of a nearby renormalized trajec- 
tory and so infer the existence of a new fixed point. One interpretation of the 
evidence for scaling is that a new renormalized trajectory exists with exponent 
K ~ 4 ~ 5, and that a new fixed point lies somewhere in the complete space of 
coupling constants with projection onto the (/3,/i) plane of ((3* ~ 7,/i* ~ —0.28). 
Evidence presented in section |5.3| shows that there is also a 'vorticity' fixed point 



associated with the 'vorticity' critical point located at about (f3 c ~ 7, \i c m —0.26) 
which terminates the first-order line. Because of the proximity of these two points 
to each other it is tempting to identify the new fixed point with the 'vorticity' 
fixed point. However, there is no direct evidence that this is so and there are 
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arguments against such an identification. The first is that we would expect the 
continuum limit defined at the 'vorticity' fixed point to be ising-like since the 
order parameter is based on a locally discrete variable: the plaquette operator, 
P x (a). The exponent k for ising-like critical points is k = 8/15 whereas for the 
new renormalized trajectory we find k ~ 4 ~ 5. This is clearly inconsistent with 
the proposal. The second argument is that the identification of the two fixed 
points means that the 'vorticity' critical point is in the domain of attraction of 
the new fixed point. The consequence is that there must be a fixed point of 
the flows in the {(A), (P)) plane in the limit that the initial lattice size is large 
enough: L —>■ oo. This is true because flows that have bare couplings held at the 
critical point values must be in the critical surface and so flow towards the new 
fixed point. In turn this implies that the correlation length £a(A/-0 for a state 
interpolated by A must diverge at the 'vorticity' critical point, (/3 c ,/i c ). This is 
unlikely since we expect £,a(/3,h) to be bounded from above by the correlation 
length in the 0(3) model at the same /3, namely £a(A °°)- This is because for 
ji < oo the presence of vortices introduces disorder in the system which acts to 
reduce the correlation length at fixed (3. However, oo) diverges only in the 

limit P — > 0, the 0(3) fixed point, and hence £a(AaO cannot diverge at (/3 c ,/z c ) 
contradicting the proposed identification of the two fixed points. Although we 
did not carry out an exhaustive investigation we found no evidence for a fixed 
point of the {(A), (P)) flows from the simulations described in section 



In figure |H] we shown an artist's impression of a possible topology of the RG 
flows in coupling constant space consistent with this interpretation and with the 
results presented in section |[ The two axes associated with the couplings /i) 
are augmented by a third which represents all other couplings. There are two 
fixed points shown. One is the usual asymptotically-free 0(3) fixed point and 
the other is the new fixed point we have identified in this work. The 'vorticity' 
fixed point is not shown. The critical surface bounds the surface of first-order 
transitions and the two phases associated with this transition are distinguished 
by the vortex density being large in one phase and small in the other. The 
line of intersection of the first-order surface with the (/3, fx) plane is the line of 
first-order transitions reported above. There are a number of possible continuum 
limits in this model, each identified with a different fixed point. A non-zero 
vortex density will be associated with the continuum limit taken at the critical 
point controlled by the 'vorticity' fixed point. At the new fixed point there are 
two relevant directions but we cannot be sure what the relevant observables are 
since in this scenario the action must be augmented by other couplings so that 
it can be tuned to lie in the critical surface and in the domain of attraction of 
this fixed point. However, the presence of this renormalized trajectory dominates 
all flows in its neighbourhood, and its influence will only be diminished if points 
in the critical surface are approached which are not in its domain of attraction. 
The example scenario of figure [10] is complicated but we have found no simpler 
toplogy consistent with the results if we demand that the scaling flows are due 
to a nearby RT in an extended model. 

A different interpretation is that the scaling flows are due to the ghost of 
the Kosterlitz-Thouless renormalized trajectory in the equivalent 0(2), or XY, 
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model. A cogent argument against a Kosterlitz-Thouless fixed point occurring in 
non-abelian models has been given by Hasenbusch in ||, but it was conjectured 
in |5| that some remnant of Kosterlitz-Thouless behaviour might nevertheless 
survive in models of this kind and give rise to the pseudo-scaling behaviour re- 
ported in II . As remarking in section [5.2| the fit to the exponent k using equation 
( |5.2.2|) should not be taken to rule out KT behaviour in favour of conventional 
second-order behaviour. Indeed, the large value for k mitigates in favour of a 
KT interpretation [14]. This explanation has the virtue of simplicity over the 
alternative picture above but it is unclear how to describe the mechanism more 
fully. 

The deviation from scaling for the large (3 points for L = 512 in sets 1 and 2 
(table |H) can be explained by noting that there is no fixed point for the ((A), (P)) 
flows and so the attempt to follow the scaling flow to larger (3 and into a fixed 
point will fail as the critical surface is approached. The conjectured renormalized 
trajectory dominates by virtue of its large exponent but scaling will eventually 
be violated as (3 increases towards (3 ~ 7. 

The strong influence of the scaling flows gives rise to a cross-over effect in the 
flows which signals the cross-over from the vortex to the spin-wave regions of the 
phase diagram. For example, in pure RP 2 this occurs at about /3 = 4.18, /i = 
0. We would naturally associate this cross-over with the observed first-order 
line but it is clear that the strength of the effect is due to the nearby scaling 
flows. This would suggest that the first-order line and the scaling flows were 
related but, as argued above, a simple relationship is ruled out and it is unclear 
whether the proximity of the two features is a coincidence or not. The region 
in which the cross-over occurs is quite narrow and has been shown as a surface 
with dotted outline in figure [H]. As (3 is increased at fixed /i through this 'cross- 
over region', the vorticity rapidly decreases from a high to low value especially 
in the neighbourhood of the critical surface. This effect means that the disorder 
also decreases rapidly and we would expect a corresponding rapid increase in the 
vector and tensor correlation lengths, £y and £t, which are deduced respectively 
from the correlators Gy(x,y) and GT(x,y) defined for RP N ^ by 

G v (x,y) = (S x -S y ) Cl 

G T (x,y) = ((S(x) ■ S(y)) 2 ) - l/N , (6.1) 

Because Gy is not gauge invariant it will vanish unless it is evaluated in a fixed 
gauge. This is analogous to the situation in QED where the electron propagator 
is not gauge invariant but the pole mass is. Technically, the gauge-fixed electron 
propagator has a cut whose discontinuity is a gauge-dependent function of a but 
whose branch point defines the gauge-invariant mass. This is due to the contin- 
uous nature of the gauge group which does not apply in our case. A reasonable 
gauge choice would be to maximize J2 x ,n a x,n ■ Gt takes the same form as the 
tensor correlator defined by Caracciolo et al. p and Sokal et al. [|J . Because 
Gt is gauge-invariant it does not require gauge-fixing before evaluation. When 
the vorticity is vanishingly small the gauge field is equivalent to a pure gauge 
and can be gauge transformed to the trivial configuration o x ^ = 1, V £E, /x. The 
physical observables in the theory are then insensitive to the chemical potential 
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\i and the theory is in the universality class of the 0(3) fixed point. 

In the 0(3) continuum limit both £y and £r will diverge, but in the continuum 
limit defined by the vorticity fixed point we expect both £y and £t to remain 
finite because, as already discussed above, the presence of disorder means that 
they will be bounded from above respectively by £y(/3, oo) and £r(/3, oo), the 
correlation lengths at the same value of (3 in the 0(3) spin model. In other 
words, at fixed (3 we expect both £y and £t to increase as \i increases, achieving 
their maximum values at \i = oo in the 0(3) model. This increase could be 
very rapid in the vicinity of the cross-over region. The operators interpolating 
the states in Gy and Gt are, respectively, Vi — Si and T t j = SiSj — l/Ndy. 
Since (Vi) = (T^-) = they show no discontinuity across the first-order line 
and hence £y and £t will not diverge at the critical point terminating the first- 
order line (gauge fixing is understood where necessary). If either of £v or £r did 
diverge it would contradict the expectation that they are bounded from above 
by their corresponding values in 0(N) as mentioned above. In principle, we 
could also study Gs(x,y) = (Uc(x)Uc(y)) c since (Uc) is continuous across 
the first-order line and it couples to the S-wave two-particle 0(N) singlet state. 
The associated correlation length £5 should coincide with £t in the continuum 
limit if the conventional scenario is assumed. We suggest that the divergent 
correlation length at the critical point is associated with the correlator of Ud(x), 
or equivalently with the vorticity correlator 

G P ((x,y)) = (P x P y ) . (6.2) 

In this study Gp was not computed. 

The pure RP 2 model (/i = 0) does not intersect any critical surface except 
the one in the basin of attraction of the 0(3) fixed point at (3 = 00 . This 
confirms the conjectures of Niedermayer et al. |7J and Hasenbusch || that RP 2 
and 0(3) have the same continuum limit. In a simulation of pure RP 2 Kunz 
and Zumbach |T5[ observe the rapid decrease in vorticity that we have associated 



with the cross-over region, and Niedermayer et al. JF| comment that in this region 
a sharp transition to a huge value for £v is to be expected. Our result is that 
the cross-over is very strongly marked in the renormalized quantities obtained 
after substantial blocking has been performed. The cross-over region separates 
two phases in one of which the vorticity density is high with a background of 
vortices pairs overlaid by a gas of free vortices, and in the other the vorticity 
density is low and does not scale as (3 — > 00. These two phases are also separated 
by a first-order line and we conjecture that a non-zero scaling limit for the vortex 
density could exist at the terminating critical point. Huang and Polonyi [TBI have 



discussed the existence of a continuum limit with a non-zero scaling vorticity in a 
generalized 2D Sine-Gordon model and the non-conservation of the kink current. 
A similar analysis could be fruitful in non-abelian models of the kind discussed in 
this paper, although it is unclear if the same techniques are directly applicable. 

In the simulation of the 2D SO (4) matrix model J| a bogus signal for scaling 
was observed which led to an incorrect measurement of the m/Ajjg ratio. In the 
context of RP 2 we would expect a similar effect for (3 ~ 3.9 because in this case 
the model renormalizes close to the scaling flows and so in the neighbourhood of 
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this coupling we should expect to see a good scaling signal. The effect is enhanced 
by the large exponent k ~ 4 associated with these flows. Simulations which are 
designed to compute m/Ajfg must have £y <C L for some largest practical lattice 
size, L. Because £y is rising rapidly in this region as a function of j3, this means 
that only a small range of f3 is useable and that this range corresponds to theories 
where the vorticity is not too low since £y would otherwise already be too large. 
The conclusion is that such simulations will see an apparent scaling due to the 
strong influence of the scaling flows. However, this scaling is not a signal for a 
continuum limit in pure RP 2 but is due to the proximity of the cross-over region 
to the scaling flows. On much larger lattices as (3 is increased a cross-over to true 
scaling would eventually be observed: the scaling associated with the 0(3) fixed 
point. However, this would be for prohibitively large values of £y, perhaps as 
large as £y ~ 10 9 0. We believe that this effect explains the results presented in 
|| who observe scaling in RP 2 (RP 3 ) but who find that the observed correlation 
length is smaller by a factor of 10 7 (10 4 ) than that deduced assuming that the 
theory is asymptotically free. We suggest that this study is actually in the cross- 
over regime where the correlation length is diminished by the disordering effect of 
vortices and the scaling, which is perhaps due to a new renormalized trajectory, 
is only apparent. The true scaling regime associated with the 0(3) (0(4)) fixed 
point will correspond to much larger correlation lengths than those studied. We 
believe that a similar effect caused the bogus signal for scaling in the analysis of 
the 5*0(4) matrix model || and the mismatch between the observed mass-gap 
and the Bethe-Ansatz prediction. 

We conclude that the RP 2 and 0(3) spin models are in the same universality 
class and that there is no evidence to the contrary. This confirms the conclu- 
sions of Hasenbusch || and Niedermayer et al. 0, but is at variance with the 
proposition of Caracciolo et al. |4j] that the continuum limits of these two mod- 
els are distinct. These latter authors propose that there is a continuous set of 
universality classes in a 2D model with mixed isovector and isotensor 0(3) spin 
interactions. The 0(3) and RP 2 theories correspond to the pure isovector and 
pure isotensor interactions respectively, and the proposition of Caracciolo et al. 
requires that these two models be in different universality classes. The work pre- 
sented in this paper shows that the opposite is true and hence that the existence 
of a continuous set of universality classes in the mixed model is unlikely. 

7 Conclusions 

In this paper we have studied the 2D RP 2 gauge model that is characterized by 
two couplings (/5,/i), where /i is the chemical potential controlling the vorticity 
computed from the gauge field plaquette expectation value. We have found that 
the role played by the vorticity in the nature of the phase diagram is crucial. Using 
standard methods we confirm the existence of a first-order transition (figures |5| 
to[7|), first suggested by Solomon et al. in the (f3,n) plane separating phases 
of high and low vorticity. The critical point terminating this first-order line 
is established to lie in the range 7.0 < f3 c < 7.5, ji c ~ —0.26, which implies 
the existence of a 'vorticity' fixed point controlling the continuous transition at 
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(PcHc). We use the Monte-Carlo renormalization group for blocking the spin- 
spin interaction and plaquette expectation values, (A) and (P), to investigate 
the topology of the renormalization group flows. We verify the presence of the 
0(3)-renormalized trajectory (at \i = oo) and find results consistent with the 
known three-loop /3-function for sufficiently large (3 once the finite lattice-spacing 
artifact has been taken into account. We establish the existence of new scaling 
flows in the ((A), (P)) plane (figures || and £|) and conjecture that they are due 
either to the ghost of the Kosterlitz Thouless renormalized trajectory in the XY 
model or to a new renormalized trajectory and its associated fixed point which 
should lie out of the (/3, fi) plane in the complete space of couplings. The scaling 
flows are consistent with a critical exponent k ~ 4 ~ 5 and the projection of 
the conjectured fixed point onto the (/3, fi) plane is deduced to be in the range 

P* « 6.5 ~ 7.5, //* « -0.31 0.26. Although the values of (&,// c ) and (/?*,//*) 

are very similar there are strong arguments against identifying the conjectured 
fixed point with the 'vorticity' fixed point. One is that the exponent k is much 
larger than that expected at the 'vorticity' fixed point, and another is that such an 
identification would imply a fixed point in the ((A), (P)) flows for bare couplings 
(/? c ,/i c ), with a consequent divergence in certain correlation lengths. This is 
contradicted by the fact that, because of the presence of non-zero vorticity, these 
correlation lengths are bounded from above by the corresponding quantities in 
the 0(3) model (/i = oo) at the same f3 which are known not to diverge for 
(3 < oo. A consequence is that the critical point at (/5 c ,/i c ) cannot be in the 
domain of attraction of the conjectured fixed point. The scaling flows dominate 
the flows in their vicinity and in particular gives rise to a cross-over (figures |6| 
and [7]) between regions of high vorticity (lower (3) to low vorticity (higher (3) 
associated by a rapid increase in the correlation length as the disorder is reduced. 
We conclude that simulations in the neighbourhood of the cross-over region for 
ji > —0.26 will show 'pseudo' scaling || because of the proximity of these scaling 
flows. The true continuum limit for such models will not be observed until true 
scaling, controlled by the 0(3) fixed point, has been established at larger (3 and 
very much larger correlation length. This is the case for the RP 2 spin model 
(fi = 0) whose continuum limit is controlled by the 0(3) fixed point and which 
is thus in the same universality class as 0(3), contradicting Caracciolo et al. 
but confirming the work of Hasenbusch [§] and Niedermayer et al. @. It 
also gives an explanation for the results discussed by Caracciolo et al. M. In 
figure [10] an artist's impression of the renormalization group flows is given for 
one scenario consistent with our results. The natures of any new fixed points 
are not established because of the known difficulty [|IJ| in distinguishing between 
fits of different scaling forms and the compatibility of the observed scaling with 
a second-order scaling form, equation Q5.2.2 ) is of phenomeno logical significance 
only. It is quite possible that any fixed point whose existence we infer from the 
data is of Kosterlitz-Thouless type. 

Our investigation has shown that the nature of gauged spin models is compli- 
cated and it is difficult to pin down more about the nature and location of the 
topological features of the renormalization group flows without more information 
concerning the relevant operators in each case. However, it is clear that a fixed 
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point in a larger coupling constant space can be close enough to the subspace 
of simple models that it very strongly influences observables and the outcome of 
tests for scaling in exactly that region accessible by simulation, namely for those 
couplings for which the correlation lengths have increased to the practical limit 
measurable on modern computers. This influence is strengthened if the expo- 
nent of the associated renormalized trajectory is large. The model studied in this 
paper is a good example of this effect. 

It would be interesting to more accurately locate the critical point at (f3 c , /x c ) 
terminating the first-order line and investigate the continuum limit it defines. 
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Figure 1: Blocking strategies for spins and gauge fields. A gauge-covariant linear 
combination of a spin and its nearest neighbours defines the blocked spin and the 
gauge field on the blocked link, which connects the solid-black sites, is assigned the 
majority sign of the three Wilson lines W+, Wo, W-. 
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Figure 2: The results of example blocking of two vortex configurations. Since vortices 
add mod2 the blocking should yield either no blocked vortices or one blocked vortex 
depending on whether the original region contained an even or odd number of vortices. 
The majority rule of equation (4.40 guarantees this important property. 
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Figure 3: The scaling flow built up from the RG flow segments of setl given in table 
||[ Each segment consists of four points corresponding to the ((A, ),{P, )) values on 
an L = 8 lattice blocked for given couplings from lattices with L = 64, 128, 256, 512, 
respectively. The segments are adjusted so that they overlay each other, and it can be 
seen from this figure, and from table ||, that the points on different segments coincide 
very well indicating that scaling holds. There are two errant points corresponding to 
the two largest (3 values and L = 512. 
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Figure 8: Two examples sets of flows, sets 3 and 4, showing cross-over of flows. Set 
3 lies to the left of set 4 and the (/?, jx) values corresponding to these sets are given 
in table Set 4 is for various in the RP 2 model (fj, = 0) in the range (3 = 3.9 to 
4.5. There is rapid variation in the renormalized ((^4, ), (P, )) values for small changes 
in (3 indicating a narrow cross-over from a region of high to low renormalized vorticity. 
For f3 ~ 4 the RP 2 flows closely follow the scaling flow of figure [3|, shown here as the 
dashed curve, implying scaling will apparently hold until, for larger /3, the cross-over 
occurs. 
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Figure 10: An artist's impression of the RG flows consistent with the simulation results. 
The 0(3) fixed point controls the continuum limit of RP 2 (fi = 0) and neighbouring 
theories. A line of critical points terminates the first-order surface and defines new 
continuum limits characterized by non-zero continuum vorticity. The observed renor- 
malized trajectory is shown associated with a new fixed point in the critical surface. 
There are indications that this fixed point does not control the second-order transition 
terminating the observed first-order line in the (/?, jj) plane, and the scenario presented 
here is consistent with this view. The cross-over region is the neighbourhood of the 
surface shown with dotted outline. 
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